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Applying a variational multiparticle-multihole configuration mixing method whose purpose is to 
include correlations beyond the mean field in a unified way without particle number and Pauli 
principle violations, we investigate pairing-like correlations in the ground states of '^^^Sn,^''^Sn and 
Sn. The same elTective nucleon-nucleon interaction namely, the DIS parameterization of the 
Gogny force is used to derive both the mean field and correlation components of nuclear wave- 
functions. Calculations are performed using an axially symetric representation. The structure of 
correlated wave-functions, their convergence with respect to the number of particle-hole excitations 
and the infiuence of correlations on single-particle level spectra and occupation probabilities are 
analyzed and compared with results obtained with the same two-body effective interaction from 
BCS, Hartree-Fock-Bogoliubov and particle number projected after variation BCS approaches. Cal- 
culations of nuclear radii and the first theoretical excited 0^ states are compared with experimental 
data. 
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I. INTRODUCTION 

Microscopic approaches based on the self-consistent 
mean field theory and its extensions are among the most 
powerful methods of describing many-body interacting 
systems. These approaches have been used for many 
years in nuclear physics [l| as well as in atomic and 
molecular physics [3,y|. In nuclear physics, they are usu- 
ally based on energy density functionals built from phc- 
nomenological parameterizations of the nucleon-nucleon 
effective interaction as the Skyrme forces [J, Q or the 
Gogny interaction [0]. 

In nuclei away from closed shells, pairing correlations 
are known to play a very important role. The techniques 
commonly used to describe them in a microscopic frame- 
work are the Hartrce-Fock (HF) -I- BCS and Hartree- 
Fock-Bogoliubov (HFB) approaches. Whereas these ap- 
proaches have proved to provide an excellent description 
of observables associated with pairing such as gaps or 
odd-even mass differences in superfluid nuclei, they suffer 
from the defect that pairing correlations are introduced 
by means of a wave-function- the BCS wave- function- 
that does not represent a definite number of nuclcons. 
As a consequence, strongly paired nuclear states contain 
spurious nucleon number fluctuations that may be large. 
In addition, part of the correlations induced by the nu- 
clear pairing interaction is not described in weakly cor- 
related regimes [ll|, [l2| . 

This problem, which arises from the particular form 
of the BCS wave-function, has led to a revival of the 
study of pairing correlations in atomic nuclei in the last 
five years. Namely, methods have been implemented by 
several groups based on particle number projected BCS 
wave- functions [ij, [ij, [I^, [la, [l3l- Whereas projection 
of the self-consistent HFB or BCS wave-functions - the 



projection after variation (PAV) technique - allows one 
to restore nucleon numbers and to compute the corre- 
sponding correction to the total binding energy, only the 
variation after projection (VAP) procedure, cither in the 
form of the approximate Lipkin-Nogami technique or us- 
ing the exact formalism is able to describe correlations 
in situations where BCS or HFB pairing is small. It is 
shown in Ref. [l^ that the two kinds of approach lead to 
significant differences in the correlation content of pro- 
jected wave-functions. 

In this work, we envisage an alternative to these pro- 
jection methods by applying a variational multiparticle- 
multihole (mp-mh) configuration mixing technique. This 
approach is similar to the MultiConfiguration Hartree- 
Fock method (MCHF) [31 well-known in atomic physics 
or to the MultiConfiguration Self-Consistent Field 
method (MCSCF) [3| employed in molecular physics, and 
it can be used to describe not only pairing correlations 
but also other kinds of long range correlations such as 
those associated with collective vibrations. The wave- 
function of the system is assumed to be a superposition 
of a finite set of Slater determinants which includes a 
HF-type state together with multiple particle-hole (ph) 
excitations of this state. Both the configuration mixing 
coefficients and the single-particle states are determined 
in a self-consistent way from a variational procedure. 

Let us emphasize that, contrary to the wave-functions 
used in the large scale shell model approach [1^, [I^ , ph 
excitations arc not restricted to those within one major 
shell. They are built from the full (finite) single-particle 
spectrum obtained in the HF-like calculation. On the 
other hand, only excited configurations involving a rela- 
tively small number of ph excitations have to be taken 
into account. This is due to the fact that single-particle 
states are derived from a self-consistent mean-field that 



already contains a large part of the effect of two-body 
interactions. It has been shown in Ref. [20[ that pairing 
correlations in usual superfluid nuclei can be accurately 
described using excitations involving the excitation of no 
more than three conjugate pairs of like nucleons. 

One important advantage of such an approach is to 
describe correlations in a formalism that explicitly pre- 
serve particle number conservation and never violates the 
Pauli principle (contrary to e.g. Random Phase Approx- 
imation (RPA) correlations). In the case of pairing, as 
will be seen in Section [Til the mp-mh method is more 
general than the fully VAP procedure. It is therefore 
in position to describe both strong pairing correlations 
without particle number violation as well as the particu- 
lar correlations occuring in weakly paired systems. 

Pioneering work along this line have been developed in 
nuclear physics using an approach referred to as Higher 
Tamm-Dancoff Approximation (HTDA) [2l|. In this 
kind of approach, the nuclear mean- field which provides 
the single-particle states is derived from an energy den- 
sity functional built with an effective force of the Skyrmc 
family, and correlations are generated by means of a 
simplified phenomenological residual interaction in the 
form of a contact two-body force. This method has been 
used to describe the behaviour of nuclei as a function of 
quadrupole deformation and the properties of long-lived 
nuclear states such as isomeric states in ^^^Hf nucleus. 
Extensions of this work, where the residual interaction is 
treated in a perturbative way, in the spirit of highly trun- 
cated shell model, can be found in Ref. [23]. An attempt 
in this direction was previously proposed in Ref. [23[ . Let 
us also mention a similar approach proposed in Ref. [24 1 



in order to describe pairing correlations in a fully particle- 
number conserving way. 

In the present mp-mh configuration mixing approach, 
the ground state and first excited states are derived vari- 
ationally from an energy density functional taken as the 
mean-value within the mp-mh wave-function of the ef- 
fective Hamiltonian built with the Gogny force. Calcula- 
tions are performed using an axially-symctric representa- 
tion of single particle states. The density distribution en- 
tering the density-dependent term of the Gogny force has 
been taken as the one built from the correlated mp-mh 
wave-function. This prescription has been adopted in the 
VAP onto particle number procedure [ij, [ij, [I^l , in the 
context of the projected HFB approach. Although there 
is no justification for using this correlated density distri- 
bution in the effective force, such a choice has been made 
in order to naturally obtain in the variational procedure 
the so-called rearrangement terms that are known to pla' 
a crucial role in the matrix elements of the mean-field 
and also of the residual interaction. The importance of 
including such rearrangement terms is clearly exposed 
and illustrated in Refs.0, [^, [^, [l3| in the context of RPA 
and QRPA methods. Let us mention that the use of a 
unique effective two-body interaction for deriving in a 
unified way both the single-particle structure associated 
with the nuclear mean-field and the correlations beyond 



the mean-field is clearly an advantage in the context of a 
completely microscopic description of nuclear states. It 
reduces the phenomenological part of the present nuclear 
structure approach as the only parameters are those of 
the nucleon-nucleon effective interaction. 

In the present work, we will only study correlations of 
pairing-type, leaving for further work other kinds of cor- 
relations. The mp-mh trial wave-function will therefore 
be restricted to a superposition of configurations taken 
as BCS-type pair excitations. Here, a BCS-type pair is 
defined as two protons or two neutrons in time-reversed 
states. By making this choice, only the usual pairing part 
of the residual interaction - the singlet even component - 
is taken into account in the description of correlations. In 
particular, proton-neutron pairing-type correlations are 
not taken into account in the present study. 

The aim of the present work is to analyze pairing-type 
correlations in nuclei in three different situations: large 
BCS pairing, medium BCS pairing and no BCS pair- 
ing by analyzing three well-known nuclei: ^^^Sn, ^"^Sn 
and ^"^"Sn. The quantities that will be examined are the 
total correlation energy, the structure of the correlated 
ground state wave-function, in particular the respective 
weights of the HF-type wave-function and of the different 
mp-mh pair excitations, and the influence of correlations 
on single-particle energies and occupation probabilities. 
Comparisons will be made with usual HFB results de- 
rived from the same Gogny two-body effective interac- 
tion and also with those of projection methods such as 
PBCS after variation. Nuclear radii and the energy and 
structure of the first theoretical excited 0"*" state will be 
compared with experimental data. 

The paper is organized as follows. In Section II 
we present the variational mp-mh configuration mixing 
method together with its restriction to the case of pair- 
ing correlations. Results obtained for pairing- type cor- 
relations in ^^^Sn, ^°^Sn and ^°°Sn are presented and 
discussed in Section HI. Summary and conclusions are 
given in Section IV. 



II. GENERAL FORMALISM 

In this part, we present the derivation of the varia- 
tional mp-mh configuration mixing method applied to 
the ground state description of even-even nuclei. We 
have considered pertinent to detail it here first as this 
approach is not commonly used in nuclear physics, con- 
trary to atomic and molecular physics. Moreover, some 
features of the method arc specific to nuclear physics, 
such as the existence of two kinds of particles and the oc- 
curence of rearrangement terms coming from the density- 
dependence of the phenomenological effective nucleon- 
nucleon interaction. 

The variational mp-mh configuration mixing method 
is a self-consistent approach that generalizes the usual 
density-dependent Hartrce-Fock (DDHF) approach [y, 
[26[ in order to take into account various types of nuclear 



correlations beyond the mean field in an unified way. 

The trial wave- function |^) that describes nuclear 
states is taken as a linear combination 



I*) = Yl ^c,„aA(t>aAau) 



(1) 



of direct products 



|(/>a„0Q„} = |</'a^)-|0a,) (2) 

of proton and neutron Slater determinants |(/)q^) and 
\4'a^). The indices tt and v stand for proton and 
neutron, respectively. Each Slater determinant |0q^), 
r=7r, V, is a multiple particle-hole (p-h) excitation a-r— 
{pihi,p2h2, ■ ■ ■)t of a HF-type reference state \(f>r) built 
with orbitals aj 



TJ- 



Ma 



i h 

where the indices h (resp. p) denotes occupied (resp. 
unoccupied) orbitals in \<J)t)- 

In Eq.ll]), the Aa^a^ are mixing coefficients. One no- 
tices that they are not taken as products of a proton 
and a neutron coefficient. The splitting of the mixing 
coefficients into the product of a proton and a neutron 
coefficient only occurs when the proton-neutron residual 
interaction is neglected. Therefore, in the most general 
case, 1^) is not the direct product of a proton and of a 
neutron wave-function. It assumes the most general form 
compatible with the separate conservation of proton and 
neutron numbers. In Eq.Q, Mq^, , indicates what we 
will call the excitation order of the Slater \(f>ar), that is 
the number of p-h excitations applied to |(/'t)- The sum- 
mation in Eq.([T]) includes the HF-type reference state 
which is obtained for M^^ = Ma^ = 0. The p-h excita- 
tions are restricted to those combinations conserving the 
quantum numbers associated to the symmetries imposed 
to the nuclear wave-function. In the present work, they 
have been taken as the parity symmetry and the axial 
symmetry around the Oz axis. Also, a finite number of 
unoccupied p- states are taken into account. Therefore, 
the number of configurations included in Eq.((T|) is finite. 

The state ([1]) depends on two sets of unknown quanti- 
ties which are taken as variational parameters: the mix- 
ing coefficients ^Q„a„ and the single particle states a^- 
entering the Slater determinants of Eq. ^ . They are de- 
termined by applying a variational principle to the energy 
functional: 

= J2 ^*o.^o.A« (0a>aJffW - A|(/)„;0„. ) (4) 

The operator H[p\ is the many-body Hamiltonian built 
with the two-body effective interaction {J12: 






^ Y{ij\v[p\\kl)ala+aiak 



ijkl 



and the term proportional to A is introduced in order to 
fix the normalisation of |^). The interaction Vij is sup- 
posed to depend on the neutron-f-proton nuclear density 
distribution p{r). Hence the notations v[p] and H[p]. As 
mentioned in the introduction, the density p(r) used in 
the two-body interaction will be taken as the one-body 
density p{r) associated with the correlated wave-function 
I*): 



p(r} = {'i'\p{rm 



(6) 



with p{r) = y 5{r — fi). This prescription is arbitrary 

since there is no justification for employing the density of 
the correlated state |^) within an energy density func- 
tional which is originally defined in the context of the 
mean-field theory. However, as pointed out in the intro- 
duction, this choice has the advantage of introducing in 
a natural way the rearrangement terms that are essential 
for obtaining realistic matrix elements for the mean-field 
and for the residual interaction 0, 0, [^, [^, [lO] . Let us 
mention that the present prescription is consistent with 
the one adopted in the application of the VAP procedure 
ofRef. [ll,[l3,[i3. 

By performing independent variations of the mix- 
ing coefficients ^a„a„ and of the single-particle wave- 
functions (flTj associated to the operators a^ , one gets 
the two extrema conditions: 




ip-rj fixed 



(7) 



= 



fixed 



Following the derivation of Appendix [^ the first condi- 
tion (O leads to the secular equation: 



'f'a.Trau 



A > , — \A 



(8) 



where the Hamiltonian matrix H, is defined by 

mnr 

(9) 



with 



and 



dV 



3?;„, = / d^r>;,„(f)^™(f) (*l^^l*) (10) 



k + v[p] 



(5) 



ijkl 

As Eq. fTU)) shows, the quantities SR^n are the matrix el- 
ements of a one-body Hamiltonian. 



One now expresses the Hamiltonian H[p\ as the sum of 
proton, neutron and proton-neutron contribution: 

H[p]=H^[p]+H''[p] + V^''[p] (12) 

Then, Eq.® takes the following form: 

mn 

E^<"^(('^<I^''MI'^".^) (13) 

+ ^KnniP,'^) {<f>a'Jatnan\(l)a^)) + 
mn 

In Eg. p^ . (7 denotes the two-body correlation function 
defined by Eq. (fT9l) below. 

From Eq. (fT3l) . one sees that the Hamiltonian matrix TL 
contains three different contributions: 

• a proton contribution involving configurations 
\4>a,4)au) ^'^d |(/)q,' (/)„' ) with the same neutron con- 
tent: |0aj = |(^a^) 

• a neutron contribution involving configurations 
\4>aT,4'a„) and |(/)a' 4>a'^) with the same proton con- 
tent: |0Q^) = I (/)„;) 

• a proton-neutron contribution. 

The first two contributions include rearrangement terms 
^mn coming from the density-dependence of the effective 
force used. Because of these terms and of the depen- 
dance of H'^[p] on the density p, the secular equation 
([5]) is a highly non-linear equation that does not reduce 
to the simple diagonalization of the Hamiltonian matrix 
TL. We give further on more details about the way this 
equation can be solved. 

From Eq. (|13p . one sees that in the variational mp-mh 
configuration mixing method, the residual interaction 
has two components. The first one originates from the 
matrix elements (0q/ \H'^[p\\<j)a^) between configurations 
\(t>a^) and \(j)ai ) differing by 2p-2h excitations and from 

the matrix elements {4'a'^4>ai\V'^''[p]\4>a^'t'a^)- The 
second one is composed of rearrangement terms between 
two configurations differing by Ip-lh excitations. 
Eq.® will be used to calculate the mixing coefficients 
^q;„q;„ for known orbitals a^ and the total energy 
associated with the correlated state. 

The second condition of ([7]) will serve to determine 
the representation used in the Slater determinants 
([3]). This condition applies because only a finite set of 
single-particle orbitals and a truncated excitation order 
are used and, therefore, the correlated wave-function 
1^) spans only a restricted part of the full many-body 



Hilbert space. Expanding the a^- over a given fixed 
single-particle basis denoted by c^: 



T] 



7 , ^n,TJ C„ ■ 



(14) 



the variation with respect to {^rj} is equivalent to the 
variation of the coefficients Cn,Tj- Hence, the second 
equation ([7]) is equivalent to: 










(15) 



As shown in Appendix |b1 this leads to the condition: 

{nH + jm^\^)md'r, atai ]|v&) =0 (16) 



where 



^*S"*^ = i ^ {mn\^\pr){^\a+a+ara,m. 

r\ / mnpr 

(17) 
Using the following definitions for the one-body density 
matrix p associated to the correlated state |^): 

p.,, = (vKla+a,!*) (18) 

and the two-body correlation matrix a: 

cTijM = (*l«i^a^a;aj|*} - p^ipik + PjkPu (19) 



one can show that Ea. (|16p is equivalent to the inhomo- 
geneous HF-type equation: 



[h[p,a],p]^Gia). 



(20) 



Here, h[p, a] is the one-body mean-field Hamiltonian 
built with the one-body density p and the two-body cor- 
relation function a: 

K, [p, a] = {i\K\:]) + r,, [p] + ary [p] + ar,, [a] (21) 

with 

T,M^Y.^My[p\\Jn) Pnm (22) 



^'^^AP\ ^ 2 ^ ("^"-1^— IW) PpmPqn (23) 



and 



m.n.nn '-^ 



Gki{a) = -Y^{im\V[p]\kn) au^r. 



(24) 



LIILIL 



(25) 



From Eq. (|25|) . one sees that Gki is an antihermitian ma- 
trix. 

The fourth term on the right hand side of Eq. ([2T|) is un- 
usual in the definition of the mean-field. It is a rearrange- 
ment field that makes the nuclear mean-field dependent 
not only on the one-body matrix p but also on the cor- 
relation matrix a. Eq. (j20p shows that the single particle 
orbitals aj' depend on a also through the matrix G{a). 
By introducing the "natural basis" associated with the 
mp-mh wave function |\1/), i.e. the representation |/i) 
that diagonalizes the one-body density matrix p: 



P^if '^/ii/^/i 



(26) 



this equation can be cast into the form: 

h[p,a]\^^)=Y,WuW^, + \XM) (27) 

where |X^(cr)) depend on the matrix G'(cr). This equation 
has the same structure as the partial differential equa- 
tions which are solved in atomic and molecular physics. 
Eq. (f27l) can be used to determine the natural states |^). 
Then, the single particle states |z) can be derived since 
they are related to the \p) by a unitary transformation 
depending only on the mixing coefficients A. 

In the first application made in this work, we have not 
solved the full equation ([77|) because of the complicated 
structure of G{a). Namely, we have neglected the de- 
pendence of h[p^ a] on a and omitted the last term in 
Eq. (|27p . With this approximation, correlations infiuence 
the single particle states \i) only through the one-body 
density matrix p. Since the simplified equation (|27p is 
still non linear, the states \p) are obtained using an iter- 
ative procedure: One starts from a HE calculation that 
gives a first set of single-particle orbitals. With the corre- 
lated wave-function |^) solution of Eq.®, one calculates 
the correlated one-body density pij ~ {^\a^ai\^) that is 
then used to calculate h. The diagonalization of h gives 
a new set of single-particle orbitals. With this new set 
of orbitals, one solved again Eq.® and the procedure is 
applied until convergence. The convergence is obtained 
when the variation of all matrix elements of p between 
two iterations is less than a defined accuracy. 

The existence of Eq. (PD|) , even when approximated by 
the scheme outlined just above, is a very important fea- 
ture of the mp-mh formalism. It expresses the fact that 
the single-particle states entering the definition of |^) 
depend on the coupling between the HE-type ground 
state and mp-mh excited configurations. Therefore, the 
self-consistent single-particle orbitals incorporate a part 
of the residual interaction beyond the usual HE single 
particle potential. This should have the consequence of 
minimizing the effects of the residual interaction and to 
allow one to truncate the expansion of \'^) to low p-h 
excitation order [20, [2l|, [23 . In this sense, the single- 
particle structure derived from Eq. ipp)) appears as the 
most adapted for describing both the mean-field and the 
correlation content of j^P). Let us note at this stage 



that the strong short range correlations due to the re- 
pulsive core of the nucleon-nucleon interaction are al- 
ready absorbed in the phenomenological effective inter- 
action. Therefore, only correlations associated with long 
range correlations have to be included, explaining why 
the above mentioned truncation in the expansion of the 
wave-function can be made. These remarks are of course 
crucial in view of the tractability of the mp-mh configu- 
ration mixing method. 

It must be emphasized that the mp-mh configuration 
mixing method does not make use of an inert core. This 
is an important feature of the method because the matrix 
elements that couple occupied deep single-particle states 
in I0T-) with high- lying unoccupied ones are not negligi- 
ble. For instance, taking the Fourier transform of the 
central part of the Gogny force shows that occupied sin- 
gle particle states can couple to unoccupied states up to 
80 MeV excitation energy. In fact these matrix elements 
are expected to contribute significantly to the renormal- 
ization of the single-particle states due to correlations. 

Since the correlated wave-function |^) is derived from 
a variational principle applied to the total energy of the 
system, it will describe the state having the lowest energy 
for a given set of quantum numbers. Then, the mp-mh 
formalism formulated here can be applied to the descrip- 
tion of ground states as well as of yrast nuclear states. 
These will be obtained as the solution with the lowest 
eigenvalue of A in Eq.®. Excited states |^i) having 
the same quantum numbers can be obtained by adding 
constraints — Ai(\E'o|^I') to the functional J- of Eq. @ 
whose purpose is to impose to |^) to be orthogonal to 
the ground state |^o) and more generally, a set of con- 
straints — X]j=o '^j(^jl^) where the |^j) are the ground 
state and excited states with lower energy than the ex- 
cited state 1^) which is looked for. This kind of extension 
will not be further discussed in this paper. A reasonable 
approximation of low energy excited states |^i) should 
however be obtained by taking the second, third, . . . lower 
energy solutions of Eq.®. One expects this approxima- 
tion to be good if the single-particle structure associated 
with 1^1) is not very different from the one associated 
with 1^0) • 

Two-body residual correlations are introduced in the 
mp-mh configuration mixing method from matrix ele- 
ments {4>ai (pa' \y[p]\4>a^4'a^) appearing in the right hand 
side of Eq.® between configurations that differ from 
two particles in two different orbitals. These matrix el- 
ements can be represented by Feynman diagrams |28| 
as in FIG. [H where the total order of excitation of 
the configuration \4>a^4>a^) {\4>a' 4>a'^)) is denoted by n 
(m). In all diagrams, p (h) stands for particle (hole) 
states. The evaluation of the many-body matrix ele- 
ments {4>a' 4>a' \y[p\\4>a^4>a„) in terms of the two-body 
matrix elements lead to three non-trivial cases for : 
a) |n— to| = 2, b) |n— 7n| = 1 and c) |n— ttiJ = 0. In the 
three cases, the two many-body configurations ](/)«„ 0q„) 
and \<f>a' 4>a' ) have to differ exactly by 2p-2h excitation 
otherwise the matrix element is zero. 
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FIG. 1: Configuration mixing diagrams 



Diagrams (al) and (a2) arc those that introduce correla- 
tions in the mp-mh configuration mixing wave-function. 
In tlie description of ground states, they mix in partic- 
ular the HF-type reference state with 2p-2h configura- 
tions. More generally, they couple mp-mh with (m-|-2)p- 
(m-(-2)h configurations. Diagrams (al) and (a2) are those 
responsible for ground state correlations in the RPA the- 
ory, where they generate virtual 2p-2h excitations (sec for 
example Ref. [i|). In the case where p2 and /12 are the 
time-reversed states of pi and hi, respectively, these dia- 
grams create Cooper pairs from the non-correlated state. 
Diagrams (bl) and (b2) arc scarcely introduced in micro- 
scopic approaches. They describe the influence of a p-h 
pair annihilation (creation) on the propagation of a hole 
(particle). They allow one to introduce the coupling be- 
tween individual and collective motion i.e., the so called 
particle-vibration coupling. 

Diagrams (cl), (c2), (c3) and (c4) that appear in the mp- 
mh configuration mixing approach between Slater deter- 
minants with the same order of excitation, are encoun- 
tered in various approaches. On the one hand, diagrams 
(cl) and (c2) are characteristic of RPA-typc correlations. 
They arc introduced in ph-RPA through the well-known 
A submatrix of the RPA matrix (see for example [l|). Di- 
agram (cl), represents the direct part and diagram (c2) 
the exchange part of the same two-body matrix element. 
Diagram (cl) describes the annihilation of a p-h pair and 
the creation of another one. In diagram (c2), a p-h pair 
is scattered from one state to another one. On the other 
hand, diagrams (c3) and (c4) appear in the pp(hh)-RPA 



and in QRPA. In such formalism, they describe pairing 
vibrations for collective states [l| . In the particular case 
where particles are in time-reversed states, they describe 
the pair diffusion mechanism of the BCS and HFB ap- 
proaches. 

Let us end this section by giving the definition of the 
correlation energy we will use. The correlation energy 
Ecorr will be taken as the difference between the total 
energy £(^) ~ (^|H|^) of the correlated system defined 
in Eq. (|A20[) of Appendix [X] and the energy of the simple 

HF method E^^ ^ (chnlHlchn): 



where 



Er. 



is such that 



£{^) - £^ 



HF 



(28) 



li^l 



= 0. 



III. PAIRING CORRELATION DESCRIPTION 
USING THE MP-MH FORMALISM 

A. Residual pairing Hamiltonian 

In this part, we apply the mp-mh configuration mixing 
formalism for the description of usual proton-proton and 
neutron-neutron pairing correlations. 
As already mentionned in the introduction, the descrip- 
tion in terms of mean-field plus residual pairing Hamilto- 
nian has played a very important role in the understand- 
ing of nuclear structure and low energy spectroscopy. 
The commonly used BCS approximation or its HFB ex- 
tension, that solve this problem in an approximate way, 
suffer from defects, as for example the non-conservation 
of particle number which is in general invocated for the 
inability of such an approach to describe weak pairing 
regimes. 

In the exact solution of a pairing Hamiltonian formulated 
by Richardson [30|, eigensolutions with seniority s=:0 and 
s=2 compete energetically. Let us note that the seniority 
quantum number gives the number of unpaired nucleons 
or twice the number of broken pairs in even particle sys- 
tems. In the language of BCS or HFB approximations, 
this means that the exact treatment produces four quasi- 
particle excitations lying lower in energy than two quasi- 
particle excitations. This behaviour is never observed in 
BCS or HFB approaches where four quasiparticlc exci- 
tations arc always far too high in energy (the first two 
quasiparticlc excitatiom corresponding to pair breaking). 
Therefore, it looks interesting to exhibit the relationship 
and differences between the PBCS solution, the mp-mh 
configuration mixing solution and the exact solution of 
the pairing Hamiltonian. 

In a standard way, when no proton-neutron pairing cor- 
relation is included, the total BCS wave-function is ex- 
pressed as the direct product of a proton by a neutron 
BCS wave- function: 



\BCS) = \BCS).^.\BCS),.. 



(29) 



In second quantization, for even-even nuclei, each BCS 
wave-function \BCS)t- (t = tt or v) is written as: 



\BCS)r^Mre^-\Q)r, 



(30) 



where |0)i- represents the proton or neutron vacuum. 
In Ea. ((30)) . Mr is the normahzation constant and B^ is 
a collective pair creation operator: 

K^W cosOrj, B+ = Y, ^aOrjbt, (31) 

with cosOrj = u-rj and sinOrj = Vrj (urj and v-rj being 
the usual variational parameters of the BCS approach). 
The 6^ operator represents the pair creation operator: 






(32) 



where the {otj, cl^j} are defined in Ea. (|14[) . In this work, 
a pair of nucleons is defined as two nucleons in time- 
reversed states and it is coupled to a total angular mo- 
mentum projection and parity K^ = 0"*". 
Expanding the exponential of Eq. ([30|) . \BCS)r can be 
written: 






N=0 



Nl 



(33) 



The \BCS)r wave- function decomposes into wave- 
functions with different numbers of particles 2N. It is al- 
ways possible to extract from Eq. ([33|) the part of \BCS)r 
having the good particle number 2N: 

0<jl<J2<---<JN 

(34) 
Defining the HF-type state with 2N particles as: 



N 



\HF)r - n ^^JO)' 



(35) 



h=l 



one finds: 



n = 0<pi<-.-<p„ -J L ^ n 

0<hi<...<h„ (36) 

n 

n kmj \bf), 

k=l 

where Af^ = 7VV H/i ^S^rh = Ilh>o sin£rh Y[p>o cos9rp . 

The wave-function \4>2n) clearly is the projection 
of \BCS) onto good particle number 2N. One sees that 
such a wave-function is a superposition of configurations 
corresponding to excitations of nucleon pairs. Eq. pS]) 
shows that the projected BCS wave-function is a subset 
of the general mp-mh wave- function of Eq.([T]), containing 
only certain types of configurations and with particular 



mixing coefficients that are products of particle and 
hole coefficients. Because of the particular form of 
mixing coefficients, \(J>2n)t is also less general than the 
exact solution of pairing Hamiltonian formulated by 
Richardson. 

In the framework of the variational mp-mh configuration 
mixing method, one considers as trial wave-function the 
reduction of ([!]) built only with the excited pair config- 
urations that are relevant for the description of pairing 
correlations with particle number conservation. Such 
kind of trial wave-function mimics the exact solution of 
pairing Hamiltonian. Then, proton and neutron Slater 
determinants are written as: 



A/a 



^«.)= Ui^tpM^Mr) 



(37) 



Here, M^^ designates the number of excited pairs in the 
configuration \4>a^)- 

As shown in Appendix [Cl without residual proton- 
neutron interaction, the mixing coefficients Aa^a^ split 
into the direct product of a proton coefficient and a neu- 
tron one. Then, the correlated wave-function takes the 
particular form: 



where 1*^-) 



Z^ 



m = \'i't).m] 



K^^,\K^)e.ndY,\Ulf = l 



(38) 



For the description of ground states of even-even nuclei, 
the proton and neutron correlated wave- functions |^^) 
and |*[J) are coupled to Rp = 0+: 

\^')o,^\n)o^.\^:)o^ (39) 

One defines the functional JF(\|/'): 

The first equation of ([7]) is equivalent to: 



dT{^') _ p dTj^') ^ ^ 



du^r 



duj: 



(41) 



One expresses the Hamiltonian H[p\ as the sum of pro- 
ton, neutron and proton-neutron contribution: 



H[p]=H^[p] + H-^[p] + V^''[p] 



(42) 



The density-dependent term of the DIS Gogny force acts 
only between proton and neutron configurations. The as- 
sociated rearrangement term is noted SH-^^ [p] . 
Following the same method as for the general formalism, 
the variational principle yields the coupled set of equa- 
tions: 

Ef^a;K'^aJi?n'^<)+C?a,a;] = [K^SuWl^ (43) 



^t/^, [((/.^JiJ-^l^^jJ+f^r^^^^.] = {\,-£^)U2^ (44) 



Equation (j43p determines proton mixing coefficients and 
equation ^^ neutron ones. The quantities Ei, and £1^'^ 
that appear in Eq. ((43|) are defined as: 



(45) 



^-= T.UZU2'SKAH''\ct>K)^ 



Y.^u2y 



>aJV""W +5^..[p]|0a,</'a„). 



(46) 
Similar expressions for £.^ and E'^'^ in Ea. (|l^ are ob- 
tained by exchanging tt and v indices in Ea. (|45p and 
Eq. (|i5)) . respectively. 

Even though no proton-neutron residual interaction is 
taken into account, the two sets of Eqs. P5|) and (P5|) 
arc not fully decoupled because of the four terms £^, £-^^ 
£^Y and E'^Y . This means that the neutron solution de- 
pends on the proton solution and conversely. 



B. Results without self-consistency 

In this part, we discuss effects of pairing correlations 
concerning the description of ^^^Sn, ^'^^Sn and ^°*'S'n 
ground states using the variational mp-mh configura- 
tion mixing approach. The same interaction is used in 
the mean-field and the residua.1 part of the Hamiltonian, 
namely the DIS Gogny force [25 • The correlated wave- 
function contains only configurations corresponding to 
pair excitations. No proton-neutron residual interaction 
is taken into account. The residual part of the Hamil- 
tonian is defined using the Wick decomposition of the 
many-body Hamiltonian H with respect to the uncorre- 
lated state |0,r0iy)- 

In this section, we focus on the effect of the mp-mh con- 
figuration mixing. Results presented in this part have 
been obtained by solving only the first equation of ([7]) 
that determines mixing coefficients: We have performed 
one HE calculation followed by one diagonalization in the 
multiconfiguration space. 

We have been interested in the convergence properties in 
multiconfiguration space, correlation energies (j28p and 
the structure of correlated wave- functions ((38l) . 
Erom a technical point of view, an eleven shell harmonic 
oscillator basis is used to expand single particle states 
(see Eg.p^) and axial symmetry is imposed. We will call 
level the twice-degenerated axially symetric state con- 
taining two time-reversed nucleon states. Because cal- 
culations are performed in even-even nuclei for which 
RP = JP = 0+ (with K the projection of the spin J 
onto the symmetry axis), the mp-mh nuclear states are 
even under the time-reversal symmetry T. Eurthermore, 
we restrict the multiconfiguration space by imposing the 
self-consistent symmetry TH2, where H2 is the reflection 



with respect to the xOz plane. Using this symmetry, 
all matrix elements and mixing coefficients can be cho- 
sen real. Let us add that the two-body center of mass 
correction term has not been included in the effective 
interaction. 



1. Convergence properties in the multiconfiguration space 

Convergence properties, according to two criteria, have 
been studied for the description of ^^^Sn, ^^^Sn and 
^'-'^Sn ground states, namely: i) the number of single- 
particle states included for the configuration mixing, as- 
suming no core (all single-particle states under the Eermi 
level have always been taken into account) ii) the trun- 
cation in the expansion of the correlated wave-function 
according to the total excitation order M = -/1/q„ + ^Ia„ 
where Ma^ defines the number of excited pairs for each 
isospin (sec Eq.([37])). 

In TABLE m one shows, for protons and neutrons, the 
number of Slater determinants restricted to M^^ = 1 or 
2 and the total size of the multiconfiguration space as- 
sociated to M < 2. Proton and neutron valence spaces 
include the entire number of single-particle levels gener- 
ated by an eleven shell harmonic oscillator basis, that 
is 286 proton and 286 neutron doubly-degenerate single- 
particle levels. 

Columns [2 — 4] give the number of proton and neu- 
tron configurations corresponding to one and two excited 
pairs. In column [5], the indicated dimension includes 
configurations as {Opair)Tr(E){Opair)^, (Opair)^(g)(lpair)y, 
(Opair)^^ (g) {2pairs)i,, {Ipair),^ ® (Ipair)^, {2pairs)T^ ® 
(Opair)^. Dimensions associated with three excited pair 
configurations are omitted as their effect is found negli- 
gible in our calculations for the three Sn isotopes. 



Nucleus 


(Ipair),! 


{lpair)„ 


{2pairs)TT 


{2pairs)^ 


Dimension 


'°°Sn 


6500 


6500 


10 101000 


10 101000 


62 705526 


''"'Sn 


6500 


7196 


10 101000 


12 434688 


69 323385 


"«5n 


6500 


8316 


10 101000 


16 698528 


80 868345 



TABLE I: Number of proton tt and neutron v Slater determi- 
nants corresponding to one excited and two excited pairs and 
total dimensions following the criteria M = Ma^ + Ma^ < 2, 
for '°°Sn,'°''Sn and ""^Sn. 286 proton and 286 neutron single 
particle levels have been considered. 

The diagonalization of the Hamiltonian matrices of Ti- 
type (see Eq.®) is accomplished using a very efficient 
technique developed for large scale shell-model calcula- 
tions [13, [231 ■ This state of the art method and the capa- 
bilities of present-day computers allow us to diagonalize 
matrices of the size presented in TABLE|T]in a reasonably 
fast way, which is a crucial point concerning the feasabil- 
ity of the present variational mp-mh configuration mixing 
method. 

The quantity for which it is obviously natural to be 



interested in is the convergence of the correlation energy 
dig). Results are shown on FIGlDand FIGElfor ^^^Sn, 
a nucleus containing large pairing correlations. 

FIG m displays the evolution of the correlation energy, 
in absolute value and expressed in MeV, as a fonction 
of the number of proton single-particle levels (proton 
orbitals are ordered by increasing energy), for a fixed 
number of neutron configurations. All neutron configu- 
rations corresponding to 286 single-particle levels have 
been included. Full circles show the correlation en- 
ergy for configurations containing only one excited pair: 
(0pa«r)7r CSi (Opair),, and {Opair),^ ® (Ipair)^. Triangles 
indicate the result obtained with in addition two ex- 
cited pairs: {Opair)^^ (E) {2pairs),j, (Ipazr)^ (g) (Ipair)^ 
and (2pairs)^ (g) (Opair)^. On each curve the first point 
on the left is calculated for 25 proton levels for which 
only the excited pair configuration occurs. This point 
gives an estimate of the neutron contribution to the to- 
tal correlation energy since in this case only neutrons 
are excited: (Opazr)^ (g) {Opair)^, {Opair),^ (^ {lpair)i, 
and {Opair)Tr {2pairs)^. For this first point, neutron 
correlations coming from one excited pair configurations 
is ~ 2.5MeV. Neutron two excited pair configurations 
bring an additional energy of ~ IMeV. As more and 
more proton single-particle levels are included, one ob- 
serves that the correlation energy saturates. Analyzing 
the correlated wave-function shows that most of the pro- 
ton correlation energy comes from the proton single par- 
ticle states close to the Fermi level. For each curve, one 
sees that a change in the slope occurs for a number of 
proton single particle levels near 50. One obtains that 
the eigensolution (black triangles) including additional 
two excited proton pair configurations and mixed two ex- 
cited pair configurations brings an energy gain of about 
~ IMeV with respect to the eigensolution containing 
only one excited pair configurations. The contributions of 
three and more excited pair configurations are not shown 
on the figure because they are small. For instance, three 
excited pair configurations contribute less than lOOkeV 
to the correlation energy. 

Similarly, FIG [3] shows the evolution of the correla- 
tion energy as a function of the number of neutron single 
particle levels for a fixed number of proton configurations 
corresponding to 286 proton single particle levels. As pre- 
viously discussed, results are shown for wave- functions 
including configurations up to one excited and two ex- 
cited pairs. The left most point on each curve, corre- 
sponding to 33 neutron single-particle levels, now gives 
an estimate of contribution from the proton to the total 
correlation energy. Adding proton two excited pair con- 
figurations is less crucial than in the case of neutron ones 
as the correlation energy gain is only ~ 200keV. Again, 
one observes a change in the slope around 50 neutron 
single-particle levels. The slope change is much sharper 
than in FIG [21 which indicates that neutron correlations 
are stronger than proton ones. One sees that the con- 
vergence of the correlation energy as a function of the 
number of neutron levels is less good than in FIGlH This 



slow convergence of the correlation energy can be under- 
stood from the magnitude of the ranges (0.7 and 1.2fm) 
of the Gaussian part of the Gogny force. When pairing 
correlations are strong, pairs involve single particle levels 
up to ~ lOOAfe^ excitation energy. Let us note that a 
similar behavior is observed in HFB calculations. 
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FIG. 2: Evolution of the correlation energy calculated with 
the mp-mh configuration mixing method, as a function of the 
proton valence space for ^^^'Sn. 
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FIG. 3: Evolution of the correlation energy calculated with 
the mp-mh configuration mixing method, as a function of the 
neutron valence space for ^^^Sn. 

Similar calculations have been done for ^'^^Sn and 
^'"^S'n. As expected, the behaviour of the correlation 
energy as a function of the proton valence space is 
similar to the one of ^^^Sn. In contrast, the convergence 
of the correlation energy with to the size of the neutron 
valence space is somewhat different for ^'^^Sn and ^'^'^Sn. 
In the case of ^^'^Sn, convergence properties according 
to the size of neutron valence space resemble the ones 
associated to proton valence space (see FIGH]), since 
^"^Sn is a doubly magic N=Z nucleus. The solution 
containing two excited pair configurations brings an 
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energy gain of about ~ SQOkeV on top of one excitatcd 
pair configurations. This clearly indicates that two 
excited pair configurations are less important for the 
description of the ^'^'^Sn ground state than in ^^^Sn. 
^'^^Sn appears as an intermediate case between ^^^Sn 
and ^^^Sn where the magnitude of neutron correlations 
is stronger than in ^'^^Sn and smaller than in ^^^Sn. The 
slope changes observed on the left of curves such as those 
in FIGH and FIGl3] are smaller than in "^S'n. Adding 
the contribution of two excited pair configurations brings 
an additional energy around 70QkeV instead of IMeV 



Nucleus 


1 T-^total 
l-t-corr 


1 T-^neutron i 
l-t-corr 1 


^oo^n 


3.67 


1.90 


'°'Sn 


4.62 


2.88 


"«Sn 


5.44 


3.74 



TABLE II; Absolute values of total correlation energy \El°r°r \ 
and neutron contribution |SSf"r''°"| for ^°°Sn, ^°^Sn and 
^^^Sn. Energies are expressed in MeV. 



The main conclusion of this study can be summarized 
as follows: i) Most correlations come from excited config- 
urations implying single particle states close to the Fermi 
level and are provided by configurations built with one 
excited pair, ii) Two excited pair configurations are es- 
sential in ^^^Sn and ^'^^Sn ground states, iii) Three ex- 
cited pair configurations can be neglected. 



2. Correlation energy 

We discuss here the magnitude and origin of the 
correlation energy obtained in the three Sn isotopes. 
Unless otherwise mentioned, all available proton and 
neutron single particle states are taken into account (286 
doubly-degenerated levels for each kind of nucleon) and 
the multiconfiguration space includes configurations up 
to two excited pairs. 

The second column of TABLE [ill displays the absolute 
value of the total correlation energy |i?corr| for ^-^^Sn, 
^'^^Sn and ^^'^Sn. One sees that ^^^Sn is the most corre- 
lated nucleus and ^^^Sn the less one: |ii^corr| = ^AAMeV 
against |i?corr| = 3.67 MeV, respectively. The third 
column gives the associated neutron contribution noted 
l^corr'^°"l- The neutron contribution is extracted from a 
calculation where neutrons are excited whereas protons 
are in the HF configuration. When one goes from ^'^"S'n 
to ^^^Sn, the neutron correlation energy increases. Let 
us note that the usual BCS or HFB approximations are 
unable to find correlations in ^'^^Sn and more generally 
when the pairing strength is small compared to the value 
of gap between the last occupied level and the first un- 
occupied level in the HF approach. 

The difference between Ecorr and i?"^™*''"" is of the 
same order of magnitude for the three Sn isotopes, about 
1.7MeT^. This indicates that correlations coming from 
protons are more or less the same, as expected. 



In TABLE IIIIl the spin-isospin two-nucleon channels 
involved for each component of the Gogny force are re- 
called (crosses) and circles indicate the channels and 
components contributing to the correlation part of the 
wave- function ([55)) . The spin-orbit contributes to the 
(S' = 1, T = 1) channel and the density-dependent 



term acts only in the mean field part. The residual in- 
teraction coming from the two gaussians arises in both 
(5 = 0, T=l) and (5 = 1, T = 1) channels [|. 

Because our method of solving the configuration mix- 
ing equations does not allow us to extract the contri- 
bution of each term to the correlation energy, we have 
studied the influence of these different terms on the cor- 
relation content of the wave-function by removing them 
selectively from the residual part of the nuclear Hamilto- 
nian. 

First, removing the Coulomb contribution from the 
residual part of the Hamiltonian changes the correlation 
energies of the second column in TABLE [U for ^'^^Sn, 
^°^Sn and "^5^ to 2.98MeV, 3.92MeV and 4.68MeV 
respectively, that is the correlation energy decreases by 
~ 70QkeV in all three isotopes. 



ST 


8=0 T=l 


S=l T=l 


S=0 T=0 


S=l T=0 


Central 


® 


® 


X 


X 


Density 








X 


Spin-Orbit 




® 






Coulomb 


<8) 


<E> 







TABLE III; Spin-isospin ST channels present in each compo- 
nent of the Gogny force (crosses). The circles indicate the 
channels and components that contribute to the residual in- 
teraction taken into account by the correlated wave-functions 
defined in Ea. ([38)l . 



Second, the different components of the nuclear resid- 
ual interaction listed in TABLE IIIIl have been succes- 
sively removed in addition to the Coulomb contribution 
mentioned above. As a result, removing all components 
except the singlet even one (S=0; T=l) leaves the cor- 
relation energy practically unchanged, and removing the 
singlet even residual interaction completely kills the nu- 
clear correlation energy. Hence the main sources of cor- 
relations are the (S=0, T=l) channel of the nuclear force 
- the one which contributes to pairing correlations in the 
HFB approach -, and the Coulomb interaction between 
protons. 
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3. Structure of correlated wave-functions 

In order to have a measure of the amount of corre- 
lations in the mp-mh wave- function |^'), we define the 
quantity T(i,j): 



Nucleus 


T(0,0) 


T(0,1) 


T(1,0) 


"«Sn 


87.21 


8.98 


3.81 


''''Sn 


88.06 


8.95 


2.98 


^o°Sn 


92.89 


3.99 


3.12 



T{i,j) = J2 1^", 



E 



IKJ'lKf (47) 



The first and second arguments of T stand for the num- 
ber of proton and neutron excited pairs inchided in |^'), 
respectively. 

TABLE lTVl displavs the values obtained for T{i,j) with 
< I < 2 and < j < 2 for the three Sn isotopes. Using 
y T{i,j) — 1, T{i,j) is expressed in percentage. One 

3=0 

observes that the HF description is approximate even in 
^°°Sn since T(0, 0) significantly differs from 100% in this 
nucleus (r(0,0) ~ 91%). The remaining 9% essentially 
come from one-pair excitation in cither the proton sector 
or the neutron one. Two-pair correlations are negligible 
in this nucleus. As expected, the two superfluid nucleus 
wave-functions contain large contributions from one-pair 
excitations in the neutron sector (~ 25%) and, to a lesser 
extent, from one-pair proton excitation (~ 2.5%) and 
two-pair excitation (~ 3.5 — 4%). It is interesting to 



Nucleus 


T(0,0) 


T(0,1) 


T(1,0) 


T(0,2) 


T(l,l) 


T(2,0) 


"«5n 


65.38 


26.04 


4.50 


2.68 


1.23 


0.17 


i°«5n 


67.44 


25.29 


3.63 


2.54 


0.99 


0.11 


i'»5n 


90.85 


5.02 


3.70 


0.16 


0.18 


0.09 



TABLE IV: Wave-function components, in percentage, for 



'Sn, 



'Sn. 



note that, if only one-pair excitations are included in the 
wave-function l^*'), significant modifications occur to the 
above numbers. This is illustrated in TABLE El The 
T(0, 0) coefficients are seen to noticeably increase and the 
T(0, 1) to strongly decrease, especially in the two super- 
fluid nuclei. This shows that although two-pair configu- 
rations have a relatively small weight in the ground state 
wave- function, their presence strongly affects the other 
components of this wave-function. Let us mention that 
three and higher order pair configurations have no influ- 
ence on the overall structure of the wave-function. Let 
us add that many one-pair coefficients \A\'^ contribute 
with small and similar magnitudes. As can be seen in 
TABLE IIVI correlations are similar in the ground state 
description of ^°^Sn and ^^^Sn and most part of them 
comes from neutron pairing. As already pointed out, the 
proton contribution is more or less unchanged from one 
isotope to the other. The difference between isotopes 
essentially comes from the neutron part. 

In FIGHl 3' schematic representation of neutron single- 
particle states pertaining to the 50-82 major shell is 



TABLE V: Components of "''Sn, ^""Sn and '""Sn wave- 
functions including only configurations with up to one excited 
pair. 
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FIG. 4: Neutron single particle levels in Sn. 



drawn. In the case of ^'^^Sn, the neutron Fermi level 
is 2^5/2 and it is completly filled in the HF approxima- 
tion. The gap between the 2^5/2 and lgr/2 subshells is 
only l.dMeV. The excitation of a neutron pair costs 
at least ~ 3.0Afey. Let note that, in BCS or HFB 
calculations with the DIS Gogny force, pairing switches 
on when the energy gap between the last occupied and 
the first unoccupied level is of the order of ~ 3.5MeI^. 
One obtains a depletion of the r(0, 0) component essen- 
tially in favor of the T(0, 1) component which is much 
larger in ^^'^Sn (~ 25.29%) than in ^°^Sn one (~ 5.02%). 
The coupling between the configurations corresponding 
to the excitation of 2^5/2 neutron pairs to the I57/2 sub- 
shell and the HF configuration is relatively strong. Those 
twelve configurations totalize ~ 15.5% of the total wave- 
function. The three configurations corresponding to ex- 
citation of 2(^5/2 neutron pairs to the 3si /2 subshell ac- 
count for ~ 1% of the total wave- function and the six 
configurations corresponding to excitation of 2d5/2 neu- 
tron pairs in the 1^3/2 subshell for ~ 3%. All the other 
configurations (~ 70 millions as indicated in TABLE |l| 
each contribute extremely small amounts. 

In ^^^Sn, the T(0,1) component is even larger than the 
^'^^Sn one. The neutron Fermi level is the completly filled 
3si/2 subshell. As can be seen on FIGlH the 1^3/2 sub- 
shell is very close to the 3si/2 one: the gap is ~ SQOkeV. 
The lowest pair excitation energy is much smaller than 
in ^'^^Sn: ~ 6Q0keV. However, the T(0, 1) component 
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of ^^^Sn is close to the ^°^Sn one. This comes from the 
fact that the larger energy gap in ^^^Sn is compensated 
by the larger degeneracy of the I57/2 state compared to 
the 3si/2 one. 

Calculations for ^^^Sn involve a total number of con- 
figurations of ^ 81 millions. As for ^'^^Sn, one can isolate 
a few configurations with strong weights: i) the two con- 
figurations corresponding to the excitation of the 3sx/2 
neutron pair into the 1^3/2 subshell account for 8.2% of 
the total wave-function ii) the six configurations corre- 
sponding to the excitation of the 3si /2 neutron pair into 
the I/111/2 subshell totalize 2.5% iii) sixty seven config- 
urations , most of them being of one excited pair type, 
have individual weights larger than 0.05% and a summed 
contribution of 12.5%. 

In the ^'^^Sn and ^^^Sn wave- functions, the two excited 
pair configurations are more important than in ^^'^Sn, 
more particularly r(0, 2) and r(l, 1) that arc equal to 
~ 2.6% and ~ 1% respectively. Comparisons of TABLE 
|TV]and TABLE El for ^"^Sn and ^^^Sn are clearly con- 
sistent with this result as the r(0, 0) and r(0, 1) com- 
ponents are strongly affected by the removing of two ex- 
cited neutron pair configurations. One sees that there 
is a strong coupling on the one hand between HF and 
one excited neutron pair configurations and, on the other 
hand between one and two excited neutron pair config- 
urations. The T(2, 0) component appears negligible for 
the description of the three Sn ground states. 

To conclude this section, let us discuss the separate 
proton and neutron contributions to the correlated wave- 
function of ^'^^Sn, the nucleus having the most correlated 
ground state. As proton-neutron pairing is not taken 
into account, the T{i,j) quantities decompose into the 
product of a proton and a neutron contribution: 



T{0,Q)=m\^.\Ulf=T^.Ti; 



(48) 



i=0,j = l 

.J^:,=o (49) 

T(l,0) = |C/;p. g \UlJ^T^,.T^ 



T(o,2) = |c/pr Y. \u:.f-n-T2 

T(2,0) = |C/;p. g \UlJ = T^.T^ 



Numerical calculations in ^^^Sn give: 



T^ ~ 93% Tf ~ 6% T^ ~ 0% 
T^ ~ 70% T'( ~ 28% T^ ~ 3% 



(51) 



One sees that, for this nucleus, even though T(0, 0) ~ 
65%, neutron mean- values of observables will be much 
more affected by correlations than proton ones. 



C. Self-consistency effect 

In this section, wc study the effect of self-consistency 
on quantities such as correlation energy, components of 
correlated wave-function, single particle spectra and sin- 
gle particle occupation probabilities. We also look at 
nuclear radii (neutron skin and charge radii) and first 0^ 
excited states, for which experimental data is available 
in most of Sn isotopes. When possible, comparisons with 
BCS or HFB approaches will be done. As explained in 
part |TT1 the full solution of mp-mh equations consists of 
solving the system of Eqs. ^ [^0]) . However, as men- 
tioned earlier (in the paragraph following Eq. (??)), in- 
stead of solving Ea. (j20p . wc have used an approximate 
procedure consisting of diagonalizing h[p\ and the contri- 
bution of a to h[p\ has been ignored. 

All the following self-consistent results have been ob- 
tained using truncated proton and neutron single parti- 
cle spaces, including the 98 lowest proton single particle 
levels and the 141 lowest neutron single particle levels. 
The total number of configurations is of the order of 10 
millions for ^^^Sn, 9 nfiUions for ^'^^Sn and 8 millions 
for ^""S'ri. One note that these numbers are significantly 
smaller than those of part IIII Bl where 286 doubly de- 
generate single particle levels were used for protons or 
neutrons. 



1. Self- consistent correlation energy 

In this section, we discuss the effect of the truncation 
of the single particle space and the role of self-consistency 
on the correlation energy. Results are presented in TA- 
BLE IVII for the four cases: mp-mh configuration mix- 
ing approach with and without self-consistency denoted 
Ef^l'; and £:^oV*r'°"S BCS and HFB approaches labelled 
respectively E^^^ and E^J^^ . Let us note that the val- 
ues indicated for BCS and HFB approximations are de- 
duced from self-consistent calculations that include the 
full single particle space associated to eleven shell har- 
monic oscillator bases (286 doubly degenerate levels). 
The BCS approximation is defined here as the reduc- 
tion of the HFB approach where only the elements of 
the pairing field matrix which are diagonal in the rep- 
resentation that diagonalizes the one-body Hamiltonian 
h[p\ are taken into account [31[. These diagonal terms 
therefore are obtained from the full Gogny interaction. 
The mp-mh correlated wave-functions include configura- 
tions built with one and two excited pairs, as discussed 
in the previous section UlI Bl By comparing E"^*^*^""* of 
TABLE Ell and £;*°*°' of TABLE HIl one observes that 
truncating the proton and neutron single particle spaces 
has a quantitative effect on the total correlation energy, 
especially in the case of ^^^Sn. The reason is that high 
energy configurations are so numerous that, even though 
they have very small individual contributions, in the end, 
they bring a non-negligible additional energy. TABLE 
IVIII shows the components of correlated wave-functions 
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Nucleus 


1 rpwith 


1 -[^without 1 


\J^corr \ 


\l^corr 


"«5n 


4.75 


3.45 


3.25 


3.86 


'"^Sn 


4.09 


3.54 


1.37 


1.73 


i««5n 


3.19 


2.79 


0.00 


0.00 



Nucleus 


T(0,0) 


T(0,1) 


T(1,0) 


T(0,2) 


T(l,l) 


T(2,0) 


"«5n 


42.09 


44.28 


3.00 


8.43 


2.09 


0.11 


i»«5n 


62.90 


28.65 


3.54 


3.62 


1.17 


0.11 


i°05n 


88.19 


6.36 


4.74 


0.27 


0.29 


0.15 



TABLE VI: Correlation energy as defined in the text for 
"^Sn, ^°^Sn and ^""Sn. The mp-mh correlated wave- 
functions including configurations with up to two excited 
pairs. Energies are expressed in MeV. 



TABLE VIIL Components of self-consistent correlated wave- 
functions for ^^^ Sn, ^'^^Sn and ^""Sn, including configura- 
tions with up to 2 pair excitation. Components are expressed 
in percentage. 



Nucleus 


T(0,0) 


T(0,1) 


T(1,0) 


T(0,2) 


T(l,l) 


T(2,0) 


"«5n 


65.06 


26.49 


4.22 


2.87 


1.21 


0.15 


i°«5n 


67.57 


26.00 


2.71 


2.84 


0.81 


0.07 


'°°Sn 


91.05 


4.98 


3.55 


0.17 


0.17 


0.08 



TABLE VIL Wave-function components, in percentage, for 
"•^Sn, ^°^5n and ^""Sn. Resuhs are deduced from trun- 
cated proton and neutron single particle spaces without self- 
consistency. 



for proton and neutron truncated single particle spaces 
without self-consistency. Comparing these values with 
those of TABLE [TVl one sees that wave- function contents 
are very similar (differences are less than 0.3%). 

From TABLE EH one sees that, comparing £'™;*^ 
and i?™^*^""**, self-consistency brings an additional en- 
ergy of the order of [400-500] keV for ^^°Sn and ^°^Sn. 
For ^^^Sn, one obtains a correlation energy larger by 
ci 1.3MeV. We will go back to this point in the next sec- 
tion. Moreover, even with smaller single particle spaces, 
the mp-mh configuration mixing approach provides sys- 
tematically more correlations than the BCS or HFB ap- 
proaches (HFB usually gives more pairing correlations 
than BCS). 

The difference between E^^^r and E^^^ is around ~ 
0.9 MeV for "^Sn, ~ 2.3 MeV for ^^^Sn and ~ 3.2 MeV 
for ^'^^Sn. The latter case is the most striking as BCS or 
HFB are not able to find correlations in this nucleus. As 
correlations coming from protons are quantitatively the 
same in the three Sn isotopes, this is a confirmation of 
the known idea that BCS or HFB approaches are good 
approximations in strong pairing regime but fail for weak 
pairing regimes. 



2. Structure of self-consistent correlated wave-functions 



We present in TABLE IVHII wave-function compo- 
nents (in percentage) obtained from self-consistent cal- 
culations. Comparing with TABLE |VlTl 7(0,0) has 
decreased by ~ 3%. In the case of ^^°Sn, this de- 
crease is counterbalanced by an increase of T(0, 1) and 
T(1,0). This is partly due to a small reduction of pro- 
ton (~ 300 keV) and neutron (~ 50 keV) gaps between 
the 2gg/2 and 2^5/2 single particle levels. The HE proton 



and neutron gaps between these two levels are respec- 
tively ~ 6.88 MeV and ~ 6.73 MeV. 

In the case of ^'^^Sn, the effect is a little more pro- 
nounced as it is accompanied by a ~ 5% reduction of 
T(0, 0). As we will discuss later, we observe a ~ 300 keV 
reduction of the proton gap between 2^9/2 and 2^5/2 and 
a ~ 70 keV reduction of the neutron gap between 2^5/2 
and 1.97/2. The HE proton and neutron gaps between the 
single particle levels mentionned just before are respec- 
tively ~ 6.21 MeV and ~ 1.86 MeV. 

In the case of ^^^Sn, the most surprising effect is the 
large depletion of the T(0, 0) component. The T(0, 1) 
component is now the largest one: c:^ 44% against ~ 42% 
for T(0, 0). We observe also an appreciable jump of the 
r(0, 2) component. As we will analyze later in more de- 
tail, this effect is essentially explained by the rearrange- 
ment of neutron single particle levels. Here again, the 
reason is the neutron reduction of gaps between 3si/2 
and 2(^3/2 and 2^3/2 and I/111/2 subshells. Consequently, 
the rearrangement of neutron single-particle states due to 
self-consistency has produced: i) a reduction of ~ 70keV 
of the 630keV gap between 3si/2 and 2^3/2 ii) a re- 
duction of ~ 300keV for the \.21MeV gap between 
2^3/2 and I/111/2 that is equal to ~ 1.2lMeV. 



Let us note that changes in single particle energies are 
bigger in ^^^Sn because the implied gaps are smaller. 

In the following, we compare mp-mh and PBCS wave- 
functions. We make use of Eg. ([55)) that gives a for- 
mal expression for PBCS wave- functions. In FIG. [5l 
the decomposition of BCS wave-functions according to 
the difference N — N' where N' is the nucleus neutron 
number is shown for ^'^^Sn and ^^^Sn. ^"^S'n, has not 
been considered since the BCS solution is identical to 
the HE solution. For ^"^S'n and ^^^Sn, the neutron 
BCS wave- functions spread essentially on eleven values 
oi N — N' . The particle number squared fluctuation 
(AiV)2 =^< BCS\N^\BCS >^ -N^ = 'iJ2k>ouM is 
6.71 for ^"^Sn and 9.95 for ^^^Sn. It is larger for ^^^Sn 
than for ^'-'^Sn because pairing correlations are stronger 
in ^^^Sn. The main component of BCS wave-functions 
in both nuclei has the good particle number N — N' = 0. 
It represents ~ 31% and ~ 25% of the total BCS wave- 
function, respectively. The components N — N' = ±2 are 
of the order of 20% and the components N — N' = ±4 
~ 10%. For ^'^^Sn, the last non negligible components 
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are the N — N' ~ ±6 ones with percentages around 2%. 
For ^^^5*71, the N — N' = ±6 components are larger and 
contribute ~ 5% whereas the N — N' ~ ±8 components 
are quite smah (~ 0.2%). 
We now look at the decomposition of the component 



Nucleus 


r(o,o) 


r(o,i) 


r(o,2) 


r(0,3) 


r(0,4) 


i°«Sn 


29.09 


54.78 


15.88 


0.25 


~0 


"«Sn 


8.03 


43.75 


38.55 


9.67 


~ 



TABLE IX: Components of ^'^"Sn and "^Sn PBCS after vari- 
ation wave-function. Components are expressed in percent- 
age. 

having the good particle number in terms of configura- 
tions characterized by a given number of excited pairs 
(see Eq. (|36p ). In order to compare this decomposition 
with the mp-mh wave-function, we have normalized this 
component to 1. Results are reported in TABLE HXl By 
comparing with the results of TABLE IVIII one sees that, 
both for ^°^Sn and ^^^Sn, the content of the PBCS after 
variation wave-function is strongly different from the one 
of the mp-mh wave-function. The r(0, 0) component is 
only ~ 29% for ^"^Sn and ~ 8% for ^^^Sn. This means 
that PBCS wave-function overestimates the contribution 
of excited configurations. 

In order to better understand these results, we have 
compared the single particle levels obtained in the three 
approaches HE, BCS and mp-mh. To this purpose, we 
have expressed the BCS single particle states \iBcs > 

h > cis linear combinations 



^mp—m, 



and the mp-mh ones 

of the HE single particle states \iHF >■ For ^'^^Sn and 
^^^5*71, we have obtained that each BCS or mp-mh single 
particle state overlaps with corresponding HE state by 



99.999%. This means that HE, BCS and mp-mh single 
particle states are very similar. The spherical single par- 
ticle states with given quantum numbers do not mix be- 
tween themselves under the effect of pairing correlations. 
Consequently, components of PBCS wave-functions de- 
pend only on the values of the variational parameters u„ 
and ii„. The difference between the T{i,j) components 
of PBCS and mp-mh wave-function appears as a conse- 
quence of the well-known fact that BCS overestimates 
pairing correlations (30| whereas mp-mh wave-functions 
are much closer to the exact ones (2011. 



3. Single particle spectra and occupation probabilities in the 
self-consistent mp-mh approach 

It is well known that, when one goes beyond the HE ap- 
proximation, the notion of single particle spectra begins 
to be lost (except in the BCS approximation). Eq. ipPl) 
illustrates this idea as h[p\ and p cannot be simultane- 
ously diagonalized. The same phenomenon occurs also 
in HEB theory [32| . However, single particle spectra can 
be obtained either by diagoanlizing the density matrix 
p and taking the mean value of h[p] or by diagonalizing 
h[p]. Here the second scheme will be used. 

In this section, we are interested in single particle level 
shifts due to pairing-type correlations. In the represen- 
tation that diagonalizes h[p] and in the special case of 
excited pairs wave-function, the one-body density ma- 
trix calculated as the mean value of the density operator 
with respect to the correlated wave-function, is diago- 
nal. Diagonal terms are directly interpreted as fractional 
occupation probabilities. 

In EIGini the neutron level of the 50-82 major shell for 
^'^^Sn and ^^^Sn deduced from four approaches, HE, mp- 
mh configuration mixing, HEB and BCS, are shown. Eor 
the two nuclei, one can note two tendencies: i) Levels 
are more compressed in the mp-mh approach than in the 
HE ones ii) The gap between the HE Eermi level and the 
next level decreases when pairing correlations are taken 
into account, the effect being the larger in HEB. 

In EIGlZl all bound proton single particle levels are 
presented in the case of ^^^Sn. ^°°Sn and ^°^Sn are very 
similar. We observe that proton single particle levels 
deduced from the mp-mh configuration mixing method 
are systematically shifted upwards. This comes from the 
well-identified effect of including the Coulomb interac- 
tion in the residual interaction responsible for correla- 
tions [33| . Besides, we obtained also a small compression 
of proton spectra as the most shifted levels are the deeper 
ones. 

In order to interpret in more detail what happens with 
single particle spectra, we examine the quantity shf — e, 
where ehf is the energy of HE single particle states and e 
the corresponding energy found in other approaches (mp- 
mh, HEB or BCS). EIG.[H] displays energy shifts of bound 
proton (upper part) and neutron (lower part) single par- 
ticle levels, between the HE and mp-mh configuration 
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mixing approaches for ^"^S'n. The vertical dashed hnes 
indicate Fermi levels. One observes that proton single 
particle states are systematically shifted upwards in the 
mp-mh configuration mixing method. One also sees that 
this shift decreases when going from the bottom to the 
top of the potential well. It is ~ 1.2MeV for the lsi/2 
state and ~ 0.5MeV at the Fermi surface. For neutrons, 
one obtains a different scenario. First, shifts are smaller 
(< 200 keV) and second shifts become positive above 
the Fermi surface. This sign inversion produces a small 
compression of the neutron spectrum. It seems intuitive 
that, at least, a part of single particle levels should be 
shifted upwards when correlations are present since the 
mean field that gives minimal total energy is the HF one. 
The level compression effects may be attributed to the 
coupling of the particle propagation with mean-field dy- 
namics, whereas the different behavior of protons and 
neutrons comes from the Coulomb residual interaction. 
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FIG. 7: ^^^Sn proton single particle spectra from the HF 
(left), mp-mh (center) and HFB (right) approaches. Only 
bounded states have been drawn. The dotted horizontal line 
is between occupied and empty single particle levels in a pure 
HF approach. 



FIG. 6: Neutron single particle spectra for ^'^^Sn (up) and 
^^^Sn (down), in the HF, mp-mh, HFB and BCS approaches. 
The dotted horizontal line is located between the HF Fermi 
level and the first empty level. 



Let us turn to ^^^Sn, a nucleus where BCS pairing 
correlations are relatively small. FIG. [5] presents the sin- 
gle particle energy shifts obtained from the mp-mh, HFB 
and BCS approaches, together with the mp-mh shifts in 
^°^Sn taken from FIG. [5] (stars). For protons, the mp- 
mh configuration mixing method predicts systematic up- 
wards shifts of the same order of magnitude or larger (for 
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the lsi/2 state) with than those of ^""S'n. The proton 
shifts olitained in ^'^^Sn by the nip-mh approach appear 
to originate from two effects: the influence of the same 
kind of correlations as in ^^'^Sn and the effect of neutron 
pairing. The differences with ^^^Sn can be explained by 
the sign of the coupling associated with pairing correla- 
tions. For example, the Isi /2 proton shifts found in the 
mp-mh approach for ^°^Sn and the HFB one for ^'^^Sn 
have the same sign, so that the total lsi/2 shift in ^'^^Sn 
is larger than the ^'^'^Sn one. On the contrary, the I/5/2 
proton shifts found in the mp-mh approach for ^'^^Sn and 
the HFB approach for ^^^Sn have opposite signs. Con- 
sequently, the final proton shift in ^^^Sn is reduced in 
comparison with the ^^'-'Sn one. 

For neutrons, a similar behavior is obtained. However, 
shifts are smaller and they appear essentially in s and p 
single particle states. 

FIG. [TO] displays energy shifts obtained in ^^^Sn. Sim- 
ilar conclusions to those of ^'^^Sn can be drawn, except 
for the 2d^/2 and l/iii/2 neutron orbitals that tend to 
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be closer to each other in the mp-mh configuration mix- 
ing approach. Besides, we have obtained for this nucleus 
an inversion between the neutron I/5/2 and 2^3/2 states 
obtained with mp-mh as well as with HFB or BCS ap- 
proaches. The HFB approximation tends to amplify this 
inversion by about ~ SOOfcey in comparison with the 
mp-mh configuration mixing method. In the mech- 
anism of the adding of the two previously mentionned 
effects {^'^^Sn type correlations and pairing), the total 
HFB or BCS shifts found for ^°^Sn and "^Sn and those 
found in ^'^'^Sn do not give exactly the shifts obtained 
in the mp-mh configuration mixing for ^'^^Sn and ^^^Sn 
but only the main trend. However, one must point out 
that the structure of orbitals is expected to change from 
one isotope to the other one. 

The influence of the residual interaction can be mea- 
sured also through the values of the single-particle states 
occupation: 

vl,=<^a+^ar.,\^>, (52) 

where |^ > represents here the ground state wave- 
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function for a given approximation (HF+BCS, HFB, mp- 
mh configuration mixing). 

In FIG Illi neutron occupation probabilities for sin- 
gle particle states located around the Fermi surface 
are drawn for ^°°S'n (top), ■^'^^Sn (center) and ^^^Sn 
(bottom). They have been calculated for five differ- 
ent cases: non self-consistent mp-mh configuration mix- 
ing (cross), self-consistent mp-mh configuration mixing 
(star), HFB (circle), PBCS (triangle) and BCS (dia- 
mond) approaches. 

In 100571, HFB, PBCS and BCS give the trivial HF 
zero or one occupation probabilities. In the mp-mh 
configuration mixing approaches (non self-consistent or 
self-consistent), neutron occupation probabilities are no 
longer equal to or 1 but they are still close to these 
values. Let us note that a similar behavior is obtained 
for protons in all three Sn isotopes calculated with the 
mp-mh configuration mixing description. Referring to 
the discussion about the structure of correlated wave- 
functions in part nil B 3[ neutron correlations can be esti- 
mated to be less than 10% of the neutron correlated wave- 
function in loo^n. Proton correlations in ^'^^Sn as well 
as in ^'^^Sn and ^^^Sn also represent less than 10%. This 
explains why occupation probabilities, in those cases, are 
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so close to the HF ones. 

In ^^^Sn and ^^^Sn, occupation probabilities for neu- 
tron are markedly different. A pronounced depletion of 
the single particle states below the Fermi level is ob- 
tained, which is compensated by a non-zero population 
of single particle states above the Fermi sea. The results 
in FIG. [H show that HFB, PBCS and BCS overesti- 
mate occupation probabilities above the Fermi sea and 
underestimate those within the Fermi sea, in compari- 
son to the mp-mh configuration mixing method. A simi- 
lar behavior has been observed when comparing HFB, 
PBCS and BCS with the results of the exactly solv- 
able model of Richardson [20, [30|- To conclude, one ob- 
serves that when pairing correlations are medium (^^^Sn) 
or weak (^°''S'n), self-consistent and non self-consistent 
mp-mh configuration mixing calculations give similar re- 
sults. In contrast, in strong pairing regime (^^^Sn), self- 
consistent mp-mh calculations give results significantly 
different from those of non self-consistent ones. Self- 
consistence induces a stronger depletion of single particle 
states inside the Fermi sea, which however is smaller than 
the one obtained with HFB. 



4- Radii 

The size and shape of nuclei strongly depend on the 
number of protons and neutrons and to a lesser extent, 
on the magnitude of correlations present in the internal 
structure. The three Sn isotopes studied in this work 
arc found spherical with the mp-mh approach. Concern- 
ing their size, wc have calculated different types of radii 
and some associated quantities directly comparable with 
experimental data in order to see the effect of pairing 
correlations obtained in the particle number conserving 
mp-mh approach: 

• the total root mean square (rms) radius: 



' T -4- T 
p ' n 



"" \^ Z + N 
with proton and neutron rms radii defined as 



(53) 



where pTi{r) and Pu{r) arc the proton and neutron 
radial densities. 

• the difference 



^'^np — ^n ^p 

a measure of the neutron-skin thickness. 
• the rms charge radius: 



(55) 



where B = Q.71AAfm comes from the proton 
form factor and 6 is a correction for center of 
mass motion. Assuming a pure harmonic oscilla- 
tor wave-function, b is given by the relation b = 
4:lA7/huj{Z+N) where the size parameter hoj is de- 
termine d by Bethe's formula hoj ~ 1.85 + 35.5(Z + 
N)-^/^ [Ij. The third contribution to re in Eg. ([56]) 
is a correction associated to neutron electromag- 
netic properties. 

The total rms radius Vav has been calculated in HF, 
HFB and mp-mh approaches for six Sn isotopes (^"'^S'n, 
^"^Sn, ^^^Sn, ^^^Sn, ^^"Sn and ^^^Sn). Results are 
shown on FIG[T21 One observes a regular increase of rav 
with the mass number A. Values obtained from HF, HFB 
and mp-mh approaches are very close to each other, the 
largest difference being found in ^^^Sn. This almost reg- 
ular increase is also obtained in the separate proton and 
neutron rms radii Vp and r„. Looking into more detail, 
one observes that for Sn isotopes containing large pairing 
correlations, Vav is larger in HF than in HFB or mp-mh. 
The same observation is true also for rp and r„ taken 
individually. This result is not very intuitive as pairing 
correlation populates levels above the Fermi level having 
on the average larger spatial extensions. However, this 
behaviour can be explained from the fact that, as shown 
previously, correlations tend to shift single particle states 
downwards, hence producing a reduction of single parti- 
cle orbital rms radii. 
One experimentally accessible quantity is the neutron- 
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FIG. 12: Total rms radii Vav for six Sn isotopes calculated 
with the HF, HFB and mp-mh approaches. Lines between 
points are drawn to guide eye. 



[rl + ^-{B'-b)-0.imi^]'/' 



(56) 



proton difference Ar„p defined in Eq. ((55l) . This differ- 
ence gives crucial indications about the distributions of 
protons and neutrons in nuclei. For that reason, it is 
a quantity more sensitive to models than the total rms 
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radius rav Results are displayed in FIG [13] for six dif- 
ferent Sn isotopes. Experimental data have been taken 
from Refs. [3J, ISa, [Sy, [33] . In the two light proton rich 
isotopes ^'^'^Sn and ^'^^Sn where no experimental data is 
available, the three theoretical approaches give negative 
values of Ar„p very close to each other. The fact that 
the proton radii arc larger than neutron ones is of course 
due to the magnitude of the Coulomb field in these nu- 
clei. In the heavier Sn isotopes, Ar„p changes sign which 
means that a neutron skin develops. One observes that 
the mp-mh approach yields values of Ar„p smaller than 
HF and HFB, especially in the most superfluid nuclei 
Ii4,ii6.i20^j^^ All calculated values are within experi- 
mental error bars, except for the mp-mh configuration 
mixing calculation of ^'^'^Sn and to a lesser extent ^^^Sn. 
The low values obtained in mid-shell Sn with mp-mh con- 
figuration mixing mainly come from the large downward 
shift of the neutron I/111/2 orbital (see FIG IlOp . which 
gives rise to a smaller value of the neutron I/111/2 orbital 
radius. One notes that the experimental error bars are 
quite large and that the different experimental results are 
scattered over a relatively large range of values. 
In order to have a more precise idea of the meaning 
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FIG. 13: Difference between neutron and proton rms radii 
Ar„p calculated with the HF, HFB and mp-mh approaches. 
Experimental measurements with error bars are also repre- 
sented. 

of results concerning Ar„p, we have evaluated also the 
rms charge radius Vc using formula (|56p . For a given nu- 
cleus, Tc depends essentially on Vp. Then, re gives an 
indication on the reliability of rp calculations. FIGfTH 
displays the evolution of re for the same Sn isotopes as 
in FIG| 13 [ Experimental data have been extracted from 
Rcfs.[3^, [3^, [40, Ell- One sees that the mp-mh config- 
uration mixing approach gives a description of charge 
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FIG. 14: Charge radii for several Sn isotopes calculated within 
HF, HFB and mp-mh approaches. Experimental data are 
indicated by red squares, with error bars. 



radii closer to experimental values, particularly in iso- 
topes where pairing is large. From this result, one infers 
that the larger deviations from experiment obtained for 
Ar„p with the mp-mh approach come mainly from the 
neutron rms radius. It is interesting to note that, al- 
though the HF and HFB Ar„p are closer to experiment, 
the separate reproduction of re and therefore of rp and 
r„ are not as good as the one obtained with the mp-mh 
approach. 



5. First excited 0^ state 

Sn isotopes manifest a very rich and complex spec- 
troscopy. In the last two decades, many experiments 
have been carried out so that, now, a lot of experimental 
data is available for most Sn isotopes between A=108 and 
A=:132, giving the opportunity to study nuclear property 
changes over a large range neutron-proton asymmetry. 
Of course, a complete description of the low energy spec- 
troscopy of these nuclei is beyond the scope of this work, 
since it would require to take into account more general 
correlations than pairing. 

Excited 0"*" states have been observed in most Sn iso- 
topes. In this study, we look first at excited 0"*" and 
compare their energy to experimental measurements. In 
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FIGfTSl the energy of the first excited O"*" state calcu- 
lated using mp-mh configuration mixing is displayed for 
i°°Sn, lO'^S'n, "'^5ri, ^^^Sn, ^^°Sn and ^^"^Sn. Energies 
vary from ~ 12MeV for the two doubly magic isotopes 
^°°Sn and "^Sn down to ~ 2MeV for mid-sheU iso- 
topes. Among the selected isotopes, experimental data 
on excited 0^ energies can be found in the literature only 
for ^^"^Sn, ^^^Sn and ^^o^^ [H. They are represented 
in FIG[T5] by horizontal red bars. The first excited 0^ 
state in ^^^Sn is known to be a collective state ( it is the 
head of a rotational band) , contrary to the second exper- 
imental excited 0^ state [4J|. Therefore, a description 
of the 0^ experimental state is beyond the scope of this 
work, and the first excited state obtained with the mp- 
mh approach is likely to be the second experimenta l On " 
excited state. This point has been stressed in Ref . [42l| . 
With this hypothesis, the difference between experimen- 
tal and theoretical excitation energies is ~ SOOkeV in the 
three isotopes where experimental data is available. One 
must note that, in the case of ^^'^Sn there is no evidence 
that the experimental first 0^ excited state is a collec- 
tive nature, as in ^^^Sn [i^]. Hence, it may be that the 
first 0+ state we calculate indeed could be interpreted 
as the first experimental 0+. For this nucleus, the dif- 
ference between experimental and theoretical excitation 
energies is ~ iOOkeV. In the case of ^^^5n, experimental 
data gives no information about the collectivity of excited 
0+ states. The experimental excitation energies of the 
three first excited 0+ states are 1.953Mey, 2.156MeT^ 
and 2A21MeV, whereas the mp-mh configuration mix- 
ing calculation give 2± 2.78MeV. 
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FIG. 15: First excited 0-1- states calculated with the mp-mh 
approach for six Sn isotopes. Experimental data for 0^ are 
indicated by horizontal bars. 



IV. SUMMARY AND CONCLUSION 



In this paper, we have presented the formalism of a 
variational Configuration Mixing self-consistent method 
adapted to nuclear structure. This approach is an ex- 
tension of the usual mean-field theory, which is able to 
treat pairing-like, RPA-like and particle-vibration corre- 
lations in a unified way. No inert core is assumed in 
this approach. In the spirit of the mean-field theory, the 
same interaction is used to describe both the mean field 
and the residual part of the effective Hamiltonian. We 
have applied this formalism to the special case of pairing- 
like correlation using the finite range density-dependent 
Gogny force. Applications to three Sn isotopes charac- 
terized by three BCS pairing regimes have been consid- 
ered: ^^'^Sn (no pairing), ^'^^Sn (weak pairing) and ^^^Sn 
(strong pairing). 

We have shown that the mp-mh configuration mixing 
method systematically finds more correlations than BCS, 
HFB and projected BCS (PBCS). In the proton sector, 
a systematic aditional energy of the order of l.TMeV 
has been found for the three Sn isotopes mentioned just 
above. For the neutron sector, we find that in strong pair- 
ing regime BCS, HFB or mpm-h provide similar correla- 
tion energy whereas in weak and medium pairing regimes 
mp-mh provide much more correlations than BCS or 
HFB. Moreover, the structure of mp-mh wave- functions 
appears quite different from those of BCS, HFB and 
PBCS. The differences manifest themselves in quantities 
such as occupation probabilities for which BCS, PBCS 
or HFB always overestimate the effect of correlations in 
comparison to the mp-mh configuration mixing method, 
and also in nuclear rms radii. 

Correlation energies have been shown to converge rea- 
sonably well using a small number of p-h excitation (up 
to 4p-4h) and a number of single particle states extending 
to ~ lOOMey above the Fermi sea. On the other hand, 
self-consistency effects, i.e. the influence of the modifica- 
tion of the single-particle states due to correlations has 
been found important when correlations are strong, e.g. 
in ^^^Sn. 

Work is in progress to include more general correla- 
tions than those considered in this study. 
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APPENDIX A: VARIATIONAL PRINCIPLE 

WITH RESPECT TO THE MIXING 

COEFFICIENTS 

In this Appendix, we derive the secular equation that 
determines the mixing coefHcients {^a„a„}- Moreover, 
we will give details concerning the evaluation of N-body 
matrix elements associated with one-body and two-body 
operators. 

The first condition given by the variational principle, ap- 
plied to the energy functional T, (see Eqs. ^ and ([7]) ) 
reads: 



a^w 

9l^ 



0. 



(Al) 



The two-body nuclear Hamiltonian is defined by: 

il = k + V\p\. (A2) 



In Ea. (jA2p . the Hamiltonian contains a kinetic term K 
(that includes the one-body center of mass correction) 
and a density-dependent potential term. As a matter of 
fact, the general formalism developed in this paper can 
be applied to any two-body interaction, as for instance 
Skyrmc or Gogny effective forces. 
Equation (jAip leads to: 



y^ Aai^c,'y\<^aAot,\'ii{p\\<\>a'^ 



E^: 



a^a^\H'a,H'a,A 



dp{r) dAl, „, 



■(fr\4>a'4>a')] 



^Aa OL,j 



(A3) 
where p[r) is the nucleon density distribution defined in 
Eq.®. 



After some manipulations, Eq. (jA3[) takes the following 
form: 



where the Hamiltonian matrix Ti is defined by 

mriT 

(A5) 
with 5R a one-body rearrangement field whose matrix el- 
ements are: 

dV 



and 



K.n = / d'r if;,„{f)iprnir) (^l^^^l^) (A6) 



5^ = iE(^^la^l^^0«.^^a. (A7) 

^ ' ijkl 



As can be seen from Ea. (jA5|) . Ti requires the evalua- 
tion of one-body and two-body matrix elements such as 
{(f>a'\o-i0.j\<t'a) and {(f>a' \cij aj aiak\(j)a) ■ In order to calcu- 
late them, excited configurations are written in the form: 



N 



i0a)=n<io) 



.10) 



(A8) 



4=1 



where TV is the number of particles (either proton or neu- 
tron) and |0) stands for the particle vacuum. The set of 
{at} indices represents orbitals that arc occupied in the 
configuration \4>a)- In order to simplify notations, pro- 
ton and neutron indices have been omitted in (jASp . The 
same notation will be used in the following. 
One assumes that, in Ea. (|A8|) . particle creation operators 
are ordered, for example by increasing single particle en- 
ergy when one goes from the left to the right. The set of 
creation and annihilation operators {af,ai} follows the 
fcrmion anti-commutation rules: 



[a. ,aj\ = aiaj+aja^ = 6.^ 



(A9) 



Using relation (|A9|) . it is easy to show that: 



N N 

m—l n—l,n^m 



Then: 



"■tajl 



N N 

U = E (-r"''5,„„a+ n <|0) (All) 

m—l n—l,n^m 



where there remains to order af within the list of aj^ 
operators. Therefore: 



N 



^+%#o) = E '^."..(-r^'+'Vr 



(A12) 



T7l = l 



with i' = i — 1 ii i < m, i' = i a i > m and l*/)^) is 
the Slater determinant obtained by removing aj from 
\(f>a), adding aj' and ordering the a"*" from left to right. 
One sees that {(j)a' \afaj\(j)a) is non zero only if |0™) and 
\(t>a') contains the same orbitals. 

For a two-body operator, the evaluation of 
{<j>a'\0'iO'to-iak\(f>a) is a little more tedious but is 
done in the same manner. One first obtains: 



iV 



aiak\<j>o.)^{l~Sik)Y.{~r+'S,^„ 



771 = 1 



N 



N 



(A13) 



E (-r'^^'^^c.^ n <io) 

with n' = n \i I < k and n' — n — 1 li I > k. 
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Then: 



a^lajaiaklf/Ja) = (1 - <^/fe)(l - S,j). 



i 3 

N 



N 






m— 1 



aja^ 



N 



n (l-'^»„)(l-<5,aja+j0) 



(A14) 



be a one-body operator. 

Using Eq. (jA12p . one obtains: 



N 



|^l|<^a)=E(E'^»')(*l^l|*) ('^"'1'^") 



TV 

EE(- 

i^j 771—1 



Z /^l 



)"+''5,o„(*|0iij)(0„,i</.r) 



(A17) 



aj a^ aiak\(j>a) = (1 - <^/fe)(l - S^j). 

N N 

m — 1 n— l.riT^rn 

(A15) 
where ( — )* and (— )■' are phases. |0J;?™") is the Slater de- 
terminant obtained by removing a J and a'^ from |(/)q,), 
adding a^^ and a'J and ordering the a+ from left to right. 
The term (0Q/|a^a^a;afc|0Q} is non zero only if If/)^™"} 
and |(/)" ) contains the same orbitals. The calculation of 
the mean-value of one-body and two-body operators is 
straightfoward using formulas Eq. ()A12p and Eq. (JAISP . 
Let 



E' 



i\9i\j)a+aj 



(A16) 



The first term on the right hand side of Eq. (|A17[) gives a 
diagonal contribution in the multiconfiguration space. It 
is a mean-field term. The second term is an off diagonal 
contribution which is non-zero only if \(l)a') = I'/"™)- In 
this case the off-diagonal term is proportional to {i\9i\j). 



Now, let 02 be a two-body operator: 



§2 = Y{ij\62\kl) afafaiUk 

ijkl 



(A18) 



Using Eq. ljAlSp . the expression of {(f)a' \()2\(f'a) contains 
three different contributions, as shown in Eq. (|Al9p . 



N 



N 



N N 

- E E(E(-)''5'^--". E ^^^cMc^'I'PT')) {^m\ki) 

N N 

' E E(E^^". E i-y-'ko^A^c'icbr)) {^m\^l) 

N N 

' E E(E(-)'"''5fc.". E hc^A4>c.'\4>T')) mo2\kj) 

N N 



N 



N 



i<j,{i,j)^k^l Kk p=l q=l,q^p 



(A19) 



The first term corresponds to the usual mean-field Let us add that the total energy £{^) of the nu- 

contribution. The four following terms as well as the last cleus is obtained by multiplying Eq. (|A4[) by A* ^ and 

one are off diagonal contributions in the multiconfigura- summing over a^ra,^- Taking into account the relation 

tion space. The four terms couple Slater determinants y^ M^ ^ p — j^^ q^iq gets: 

|(/)q) and \4>a') differing from one particle in one orbital a^a^ 

and the last term couples two Slater determinants that 

differ from two particles in two different orbitals. £:(^) = A - ^ K;,„(*|a+„a^„|*) (A20) 
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APPENDIX B: VARIATIONAL PRINCIPLE 

WITH RESPECT TO THE SINGLE PARTICLE 

ORBITALS 



We detail here the derivation of Eq. (pOl) . The starting 
point is the second condition of system ((T]) where one 
assumes fixed mixing coefficients. 






0. 



(Bl) 



Using Eq.(|3]), the variation 5T{'i') of the energy func- 
tional is equal to: 



+ {'^\SV[p]\^) 



with 



and 



SV[p] 



^3p dV[p] 
dp[r) 



5p{r) 



(B2) 



(B3) 



5p{r) = ((5*|p(f)|*) + (*|p(f)|5*) (B4) 

First, let us note that the variation of \^) with respect 
to the orbitals a^ can be written: 

\6^) = iS\^) (B5) 

where S is an infinitesimal herniitian one-body operator. 

S = Y,Suatai (B6) 

kl 

In fact, using Thouless' theorem, a variation of the or- 
bitals can be written: 



e'^a+e-'^ 



a+ + [iS, a+] 



(B7) 



Therefore, any Slater determinant of the form (jASp varies 
according to 



10a) 



e"|0«)^(l + i5)|0„) 



(B8) 



as e |0) = |0). Consequently 



|vI/) = ^A„|0„) ^ |vI/) + |5vi/) = ^A„(l + i^)|0„) 

a Q 

(B9) 



which yields Eq. fBS)) . 

Let us mention that \5'^) represents a genuine variation 

of 1^) only if |^) belongs to a subspace of the full N- 

particle Hilbert space. This is the case here since the |^) 

is built from a finite set of mp-mh excitations. 

Using (jB5p . (jB2p can be expressed as: 



5T{^) = i(*|(iJ - A + / rfV (*|^^^|*)/5(f))5' 



p{r) 



(BIO) 



The second condition ([7]) finally leads to 

,dV[p] 



mH- 



dV(*| 



9p(r 



|*)p(f), a+ai]\^)=0 (Bll) 



Let cr be the two-body correlation matrix defined as : 

(^la+a+ana/l*) = piiPnm ^ PlmPni + CTu^mn (B12) 

Eq. ljBlip can be seen to be equivalent to the equation 

[h[p,a],p]^G{a) (B13) 

with 

Gki{<j) = -^{im\V[p]\kn)aa^mn 



-^{■ml\V[p]\ni)akz. 



(B14) 



Eq. (jB13p appears as an inhomogeneous HP equation, the 
right hand side G{a) being an antisymmetric matrix de- 
pending only on the two-body correlation matrix a. This 
equation reduces to the usual HF condition when a is 
taken to be zero. 



APPENDIX C: PROTON-NEUTRON SPLITTING 
OF THE MIXING COEFFICIENTS 

In this Appendix, one assumes that the correlated 
wave- function j^P) is particular in such a way that the 
residual proton-neutron interaction part of H deduced 
from Wick's theorem gives no contribution. This is the 
case for the correlated wave-function used in the pairing 
application of part [1111 . Then, let us defined the restricted 
Hamiltonian Tirestr. containing only the terms that give a 
contribution with respect to the correlated wave-function 
that has been chosen. It can always be written as the sum 
of a proton and a neutron contribution: 



fire 



T-C + T-C 



(CI) 



In this case, Eq.® is equivalent to: 

^(0ajH"|0<) ^a;a,.+ 

^(0aj7T:n0<) ^a„< - E A, 



(C2) 



Since the matrices associated with T-f^ and T-C are her- 
mitians, they can be diagonalizcd using unitary matrices 
C/'^ and U": 






f^a^a'^^a'^^k " ^a^,k^k 



(C3) 



(C4) 
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A consequence of Eq. (|C3|) is: 
Applying Y.C. K:,j to Eq.jCl gives: 



(C5) 



(C6) 



By comparing Eq. ()C6|) with Eq. ()C4|) . one sees that, if the 
mixing coefficients A and the total energy E are solutions 
of Eq.dCH, then the quantities J^a^alj^Ka'^^ should 
be proportional to one of the UJ^, j. : 



7 , ^al,j^a'^al ~ ^jkU^'^ 



(C7) 



Using Ea. ljCTD . 



An 



J 



JJf & Tf^ y^TT TJTT 



(C9) 



Therefore A^i a' , solution of Eq. (|C2p with the eigenvalue 



Ej, is proportional to U^, f^ 



E^El 

In the same way, by exchanging proton and neutron in- 

= C^, U^, J is an eigenvec- 



dices, one shows that A^ 



tor of dell) with the eigenvalue E = E^ + E^. Then, C^^ 
is necessarily proportional to C/^, ■ . 
Taking into account the condition ^^ l^a^a^P = 1, one 
finally obtains: 



Ar, 



TJTT TJ-U 



(CIO) 



where ^jk is a complex phase, and (|C6p shows that 



e = ej: + sj 



(C8) 



Consequently, with the form (jCip Tirestr. , the mixing co- 
efficients Aa„a„ are products of a neutron and a proton 
contribution. 
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